{
 "metadata": {
  "name": ""
 },
 "nbformat": 3,
 "nbformat_minor": 0,
 "worksheets": [
  {
   "cells": [
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "%run Gas.py"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 115
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "gasHO=Gas(9)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 116
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "gasHO.vec_r[0][0]=-9.1\n",
      "gasHO.calculateDeltasHO()\n",
      "print gasHO.deltaX\n",
      "gasHO.vec_r"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "[[ 0.   0.1  0.   0.   0.   0.   0.   0.   0. ]\n",
        " [ 0.1  0.   0.   0.   0.   0.   0.   0.   0. ]\n",
        " [ 0.   0.   0.   0.   0.   0.   0.   0.   0. ]\n",
        " [ 0.   0.   0.   0.   0.   0.   0.   0.   0. ]\n",
        " [ 0.   0.   0.   0.   0.   0.   0.   0.   0. ]\n",
        " [ 0.   0.   0.   0.   0.   0.   0.   0.   0. ]\n",
        " [ 0.   0.   0.   0.   0.   0.   0.   0.   0. ]\n",
        " [ 0.   0.   0.   0.   0.   0.   0.   0.   0. ]\n",
        " [ 0.   0.   0.   0.   0.   0.   0.   0.   0. ]]\n"
       ]
      },
      {
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 120,
       "text": [
        "array([[-9.1, -0.5,  0. ],\n",
        "       [-8. , -0.5,  0. ],\n",
        "       [-7. , -0.5,  0. ],\n",
        "       [-6. , -0.5,  0. ],\n",
        "       [-5. , -0.5,  0. ],\n",
        "       [-4. , -0.5,  0. ],\n",
        "       [-3. , -0.5,  0. ],\n",
        "       [-2. , -0.5,  0. ],\n",
        "       [-1. , -0.5,  0. ]])"
       ]
      }
     ],
     "prompt_number": 120
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "gasHO.deltaR2;"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 109
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "gasHO.calculateForcesHO()\n",
      "print gasHO.Fx\n",
      "gasHO.force"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "[[ 0.   0.1  0.   0.   0.   0.   0.   0.   0. ]\n",
        " [-0.1  0.  -0.   0.   0.   0.   0.   0.   0. ]\n",
        " [ 0.   0.   0.  -0.   0.   0.   0.   0.   0. ]\n",
        " [ 0.   0.   0.   0.  -0.   0.   0.   0.   0. ]\n",
        " [ 0.   0.   0.   0.   0.  -0.   0.   0.   0. ]\n",
        " [ 0.   0.   0.   0.   0.   0.  -0.   0.   0. ]\n",
        " [ 0.   0.   0.   0.   0.   0.   0.  -0.   0. ]\n",
        " [ 0.   0.   0.   0.   0.   0.   0.   0.  -0. ]\n",
        " [ 0.   0.   0.   0.   0.   0.   0.   0.   0. ]]\n"
       ]
      },
      {
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 110,
       "text": [
        "array([[ 0.1,  0. ,  0. ],\n",
        "       [-0.1,  0. ,  0. ],\n",
        "       [ 0. ,  0. ,  0. ],\n",
        "       [ 0. ,  0. ,  0. ],\n",
        "       [ 0. ,  0. ,  0. ],\n",
        "       [ 0. ,  0. ,  0. ],\n",
        "       [ 0. ,  0. ,  0. ],\n",
        "       [ 0. ,  0. ,  0. ],\n",
        "       [ 0. ,  0. ,  0. ]])"
       ]
      }
     ],
     "prompt_number": 110
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "gasHO.calculateDeltas()\n",
      "gasHO.deltaX"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 111,
       "text": [
        "array([[ 0. , -1.1, -2.1, -3.1, -4.1, -5.1, -6.1, -7.1, -8.1],\n",
        "       [ 1.1,  0. , -1. , -2. , -3. , -4. , -5. , -6. , -7. ],\n",
        "       [ 2.1,  1. ,  0. , -1. , -2. , -3. , -4. , -5. , -6. ],\n",
        "       [ 3.1,  2. ,  1. ,  0. , -1. , -2. , -3. , -4. , -5. ],\n",
        "       [ 4.1,  3. ,  2. ,  1. ,  0. , -1. , -2. , -3. , -4. ],\n",
        "       [ 5.1,  4. ,  3. ,  2. ,  1. ,  0. , -1. , -2. , -3. ],\n",
        "       [ 6.1,  5. ,  4. ,  3. ,  2. ,  1. ,  0. , -1. , -2. ],\n",
        "       [ 7.1,  6. ,  5. ,  4. ,  3. ,  2. ,  1. ,  0. , -1. ],\n",
        "       [ 8.1,  7. ,  6. ,  5. ,  4. ,  3. ,  2. ,  1. ,  0. ]])"
       ]
      }
     ],
     "prompt_number": 111
    }
   ],
   "metadata": {}
  }
 ]
}